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Abstract 

We develop some properties of the non-BPS attractive STU black hole. Our principle 
result is the construction of exact solutions for the moduli, the metric and the vectors 
in terms of appropriate harmonic functions. In addition, we find a spherically-symmetric 
attractor carrying p° (D6 brane) and q a (D2 brane) charges by solving the non-BPS 
attractor equation (which we present in a particularly compact form) and by minimizing 
an effective black hole potential. Finally, we make an argument for the existence of multi- 
center attr actors and conjecture that if such solutions exist they may provide a resolution 
to the existence of apparently unstable non-BPS "attractors." 
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1 Introduction 



The black hole attractor story (which began with BPS states - [JJ, [2J, [3]) has experienced a 
recent flurry of activity with regards to the existence of non-BPS attractors ([J], [5]). Other 
interesting recent developments include the development of the attractor mechanism for flux 
vacua ([6], [7]), the features of Af = 8 attractors ([8], [9]) and several additional properties of 
both BPS and non-BPS attractors ([UJ, P2], [IS], [Hi], pi] and [15]). 

In particular, it seems that many of the interesting properties and features of attractive 
BPS configurations are shared by non-BPS ones, so long as the black hole in question remains 
extremal. For example the non-BPS attractor equation was proposed in [B] and developed in 

ng. 

Such recent developments lead one to speculate on what other attractive features may apply 
to non-BPS extremal black holes. However, as we try to explore this situation, an immediate 
problem arises. The first order BPS equations afford considerable simplifications that are not 
present for the more general non-BPS situation. As a result it is considerably more difficult to 
make progress understanding the properties of non-BPS attractors. One way to sidestep this 
calculational intractability is to consider a simple sub-class: The STU black holes. As outlined 
in more detail below, STU black holes can be constructed from compactifications of type IIA 
string theory whose moduli space is described by four homogeneous (or three inhomogeneous) 
co-ordinates. In particular, our goal in this paper is to demonstrate that, for the STU black 
hole, we can find exact expressions for the moduli fields, the metric and the vectors throughout 
the spacetime. 

It is well established that, along with giving the attractive values of moduli, the BPS 
attractor equation can also be used to give an exact expression for the moduli fields everywhere, 
by taking the attractive values and replacing charges with appropriate harmonic functions 
(details can be found in [TT] , [IB], [JJJ] and [2U]). While proving the equivalent statement for 
non-BPS black holes is difficult, it can checked explicitly in the STU case, where we will discover 
that the exact expressions for the moduli and the metric are given by a similar prescription as 
for the BPS attractors. 

This paper is organized as follows. In section 2, we review the non-BPS attractor equation 
for STU black holes and find a new non-SUSY solution for a black hole carrying charges 
T = (0,q a ,p°,0) (corresponding to wrapped D2 and D6 branes, and manifestly dual to the 
already known DO-DA system). Next, in section 3, we perform the calculation advertised above, 
and confirm that the non-BPS attractor has its moduli described by appropriate combinations 
of harmonic functions. Finally, we speculate on the existence of multi-centered non-BPS black 
holes. We argue that not only are such solutions plausible, but also that they will exhibit a 
similar split attractor flow to the analogous BPS situation and also allow us to resolve apparently 
unstable "attractors" (for example, non-BPS STU extremal black holes with charges T = 
{qo,Qa,P°,P a ), which have fixed values of the moduli at the horizon, but are unstable to small 
perturbations). 
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2 STU Black Holes and the Attractor Equation 



We present here a brief review of the STU black holes and attractor behavior - further details 
can be found in @], 0, [E] and [21]. 

When type IIA string theory is compactified on a T 6 (or, equivalently, when M-theory is 
compactified on T 7 ) one recovers M = 8 supergravity in d=4 with 28 vectors and 70 scalars in 
the coset space ^ = , with an on-shell U-duality symmetry £7(7). The STU model is a 
consistent H = 2 truncation of a A/" = 8 model. The model depends on three complex moduli 
(parameterizing the coset space (5L(2,R)[/(l)) d ) and has 4 vectors (one as an M = 2 partner 
of the graviton and the others in vector multiplets). 

The most general one-center BPS black hole solution of an Af = 8 model preserving 1/8 of 
the available supersymmetry is a solution of the STU model, modulo U-duality [22]. To generate 
the most general solution one has to act with H = SU(8) and transform from the 4 + 4 (electric 
and magnetic) charges of the STU model to the 28 + 28 charges of the M = 8 model. This 
corresponds to the decomposition of the 56 of £7(7) with respect to 5*0(4,4) x (5L(2,Z)) 3 : 

56 -> (8,, 2, 1, 1) © (8 S , 1, 2, 1) © (8y, 1, 1, 2) © (1, 2, 2, 2) . (2.1) 

The charges of the STU model correspond to a singlet in 5*0(4, 4) and transform as a (2, 2, 2) 
irreducible representation of (SL(2, Z)) 3 - shown in the last term of fl2.ll) . Acting with G = £7(7) 
one can generate the generic values of scalars at infinity. 

The bosonic part of the action of the four dimensional M = 2 supergravity (without hyper- 
multiplets) (in units with G^ — 1) is described by 

S= lLJ dAx ^{ - f + G «» dza ■ d ^ + Im (AA S )^ A • ^ + Re(A" AS )^ A • (*^) s ) . (2.2) 

The fermionic part (the gravitino and the chiral gaugino) can be constructed via the supersym- 
metry transformations: 

5^ = D^a + e AB T-y»e B , (2.3) 
S\ aA = z 7 ■ dz a e A + % -e AB T~ a ■ 7 e B . (2.4) 

€a is the fermionic parameter of the transformation and Eab is the 50(2) Ricci tensor. In 
above transformation laws, T~ and T~ a are the two-form graviphoton and vetor multiplet field 
strengths respectively. The integrals of these two two-forms on some two cycles give us the 
central charge of the supersymmetry algebra and its moduli space covariant derivative: 

Z = ~f T-,D & Z = ~f T~ h G h - a . (2.5) 
2 Js 2 1 Js 2 

The STU models correspond to the subset of the above theories with 3 vector multiplets. 
The moduli space of the theory is described by the four homogenous coordinates X A . These 
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coordinates combine to give the STU prepotential 0: 

yl v2 v3 

F(X) = . (2.6) 

The Kahler potential and superpotential of the theory can be readily constructed from ( 12. 6j) . 
Working with the inhomogeneous coordinates z A = ^ = (1, (o £ {1)2,3}) the Kahler 
potential is: 

K = -ln(-z(z 1 -z 1 )(^-z 2 )(/-z 3 )) . (2.7) 
The metric and connection on the moduli space then follow immediately: 



$ab ^ab _ zab(„a -a\2 — " 2 

~a\2 



Gab = — ( ~ ™ , G ab = -6 ab (z a - z a y , = - (no sum on a) . (2.8) 



Let's initially assume a static, spherically-symmetric spacetime, with a metric of the form: 

ds 2 = -e 2U{r) dt 2 + e- 2C/(r) dxW (2.9) 

In this background it is straightforward to define electric and magnetic charges (details can be 
found in appendix A of [16]) which give rise to a superpotential of the form 

W = q A X A -p A d A F , (2.10) 

where q A and p A are electric and magnetic charges respectively. The set of charges (go, Q a ,P°,P a ) 
corresponds to the charges of (-DO, D2, D6, DA) branes wrapped on (0, 2, 6, 4) cycles. 

An effective potential for the moduli (in this spherically-symmetric, static spacetime )is 
given by: 

V BH = \DZ\ 2 + \Z\ 2 = G a ~ b (D a Z)D h Z + ZZ . (2.11) 

Z the central charge of the supersymmetry algebra and, in our case, is given by Z = e K ' 2 W ' . If 
the black hole in this solution is extremal (i.e. has zero temperature) then it will be an attractor 
for the moduli - the values of the z a will be fixed at the horizon, independently of their values 
at infinity ([1], [IS]). 

The horizon values of the moduli can be obtained by minimizing the above effective poten- 
tial - either directly or through the attractor equation described below. Note, however, that 
the extremal point of Vbh must be a minimum, else the attractor will be unstable to, say, 
perturbations away from staticity. For BPS attractors this condition is always satisfied, for 
non-BPS ones, however, we have to check the second derivatives of Vbh explicitly. 



2.1 The Attractor Equation 

In this section, we derive an identity using the symplectic properties of M = 2 supergravity. 
Then, by imposing the minimization condition of the black hole potential on the identity we 
construct an algebraic relationship between the charges and the attractive values of the moduli. 

1 When string theory is compactified on some spaces which break supersymmetry down to N = 2, e. g. 
heterotic theory on K3 x T 2 , the prepotential (|2.6[) has corrections. However, when considered as a consistent 
truncation of Af = 8 model, the prepotential is given just by eq. (|2.6[) . 
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If we define the covariantly holomorphic period vector II = (L A Ma), then the symplectic 
structure of Af = 2 supergravity relates the upper and lower components of the period vector 
via the vector coupling in the following way: 

M A = A/" AS L E , D a M A = A/aeA^ , (2.12) 

where IT = (L A M A ) = e~ (X A F A ), and (X A Fa) is the holomorphic period vector. Now, we 
form a matrix by constructing the tensor product of the period vector covariant derivative with 
its complex conjugate 



G ab D a U ® DuU = G ab 



D a L A D b l* D a M A D- b L s 
D a L A D b M^ D a M A D- b Mx 



yab i D a L A D b L s M AA D a L A D b L* 



V U^D a L A D b L A M AA KT^D a L A D b L T ) ' [ZA6) 
where in the second line we have used (I2.12p . If we use the following special geometry identityE] 

G ab D a L A D b ZP = -hm(M-Y^ ~ L A Lp , (2.14) 
then we find another useful identity 

^ ~ * ( 1 A 1 



G M D a U®D l Yl=--l _ 1 q j--M-U®U, (2.15) 

where matrix M. is defined as 

/ Im(A/--r s Re (TV) a a Im ( A/" - 1 ) As \ 

71 " V Im (7V~ 1 ) A A Re (TV) ae Im(Ar) AS + Re(Ar)AAlm(AT- 1 ) Ar Re(Ar)rs J ' 1 } 

In fact, (12.151) expresses the tensor product of the covariant derivative of the period vector 
with its complex conjugate in terms of the tensor product of the period vector itself with its 
conjugate via the vector couplings. 

Now, assume that Y is the set of magnetic and electric charges Y = (p A q A ). We define Y 
by a symplectic rotation: 

_°! J ) r ' ( _°i l) eS P ( 2 ( /i1 * 1 + X )' Z ) • ( 2 - 17 ) 

Recalling (12.101) . it can easily be seen that the central charge can be expressed as Z — T* • II. 
Using this form, we can compute G ab (D a Z)D b lI using (12.151) : 

G ab (D a Z)D b YL = f ' • (G ab D a Yl <g> D b Yl^ . (2.18) 

Thus, we have the following expression: 

2iZYl + 2iG ab (D a Z)D b Yl = Y - iY f ■ M . (2.19) 



This identity can be proved by considering the inner product (D a Il, D^H) ~ ~iG a i and using (|2.12[) . 
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It has been shown in [23], and [23], that the symplectic invariant I\ is given by: 

h = \Z\ 2 + \DZ\ 2 = -i f* -M -f . (2.20) 

Finally, we differentiate the above expression with respect to the charges and substitute for 
(EI9D: 

T + i^- = 2iZU + 2iG ab (D a Z)D- b U . (2.21) 
dT 

So far, everything we have said is generically true and independent of the detailed model. 
However, if we restrict ourselves to STU models, then the symplectic invariant l\ is a function 
of charges. It is given by l\ = ^/|W(r)|, in which W is given by: 

W(r) = 4((p 1 g 1 )(p 2 g 2 ) + (Ai)(p 3 ? 3 ) + (P 2 q2)(p 3 q 3 )) - GAa) 2 -4Ai<?2<? 3 +%W • (2.22) 

In [T7], Berhrndt et al. demonstrated that W(r) is an [SX(2,Z)] 3 invariant which uniquely 
determines the form of the STU metric and further — W has been established as the Cayley 
hyperdeterminant by Duff in [25] . The fact we must use the absolute value of W in order that 
this expression be valid for non-BPS attractors was discussed in [Tj3] - here it was observed that 
W(r) is positive for BPS attractors and negative for non-BPS ones. 

Recalling that at an attractor point the effective black hole potential will be minimized (see 
[23J, [3], [IS]) we note that extremization of (12. lip gives: 

2ZD a Z + %C^rr{T)f){T) t Z) = . (2.23) 

Hatted indices are for the tangent space and rj ab is a flat Euclidean metric. Here G aa = e??7 aa e2. 
The rank 3 completely symmetric tensor C a g £ in the tangent space is related to the covariantly 
holomorphic symmetric curved space symmetric tensor C a b c as follows: C^g- = e a h ^e%C a \ )C . Here 
the rank 3 completely symmetric covariantly holomorphic tensor of the special Kahler geometry, 
C a bc, is defined as follows [2 



C abc = e K {d a X A d b X B d c X c )d A d B d c F . (2.24) 

For the STU model with the prepotential (12.61) the non-vanishing components of the symmetric 
3-tensor are Cj^ = 1 

Note that the covariantly holomorphic symmetric 3-tensor C a bc has a chiral weight 2 under 
Kahler transformations whereas in the tangent space has a chiral weight 0. 

The chiral U(l) weights of various objects living on Kahler manifolds are defined as follows. 
The Kahler potential transforms as: 

K^K + f(z)+f(z) . (2.25) 

The covariant derivative of an object with the chiral weight p is given by = (d + ipQ)§ 
where the Kahler connection is given by Q = —^(d a Kdz a — d a Kdz a ). In components this 
means: 

D a ®= (d a + ^pd a K)<f> . (2.26) 
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Thus the chiral weight of Z is 1, Z is -1, e~ is -2/3 and ei is 2/3. This means, in particular, that 
under Kahler transformations the central charge Z as well as its covariant derivative transform 
by a phase: 

Z^Ze 1 ^-, D~ a Z -f D-aZe 1 ^- . (2.27) 

Further, observe that (12.231) has a chiral weight -2/3 whereas (12.211) has a vanishing chiral 
weight. In [9] (12.231) was studied by embedding the Af = 2 supergravity into an Af = 8 one. 
Following the same method we define Y = iZ, and Y a = DgZ. This gives us the solution: 

Y Q = pe i{n - 3 ^ , Y a = pe^ . (2.28) 

Here the arbitrary phase <fi reflects the possibility to perform a Kahler transformations in 
agreement with the weights shown in (12. 27ft . (12.211) can be rewritten in the following form: 

dl 3 

r + ^ = 2p(n + ^A i ri) (2.29) 

a=l 

The arbitrary phase <fi drops from this equation since the left hand side of eq. (I2.21l is manifestly 
Kahler invariant and in the right hand side the phase cancels in each term separately. 

Now, if we divide (12.291) by its zeroth component, then we follow [9] and obtain: 

^ = + ^f 1 (2-30) 

= M A + j EL D&M A 

In fact, these equations are dual to each other and either set is sufficient to completely determine 
the moduli at the horizon. 

We now illustrate that an even simpler form of these equations can be found, as the r.h.s 
simplifies significantly. Consider (I2.30p . since (with the vielbein e? given by e? = i5^(z a — z a ))\ 

D h L c = e\D a L c = e K ' 2 e a b {5 c a + (d a K)z c ) , (2.32) 

we can readily establish: 

L c + ij^ D&L C = -2e K ' 2 z c , L° + iJ2 D&L = -2e K ' 2 . (2.33) 

a=l a=l 

As expected the above is independent of charges, and the values of the moduli at the horizon 
for any T are thus given by: 

^) = .a/% • ( 2 - 34 ) 

P — 1 a 
r dqo 

This appears to be the simplest form of the non-BPS attractor equation for the STU model. 
In order to obtain eq. (2.34), we have used the gauge in which X° was real. However, the 
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attractor equation can be generalized^] to the case in which X° = e 3 * , with real 5. In this 
gauge the attractor equation generalizing the one in eq. (2.34) becomes 



z c - e i5 (2z c + z c ) _p c + i 8?( 



aim 

(2.35) 



P + 1 dq 

and similar equation for component M\. In case that 5 = 2nir eq. (I2.35I) reduces to eq. (2.34). 
This gives a correct solution for (D2-D6) and (D0-D4) systems. However, if one considers 
(D0-D6) system, then eq. f!2.35|) requires 5 = (2n + 1)tc. 

For the BPS case Y a vanishes and so a similar procedure to that used above will result in 
an attractor equation of the form of (I2.29p . but with the r.h.s. containing only the first term. 
Dividing such an expression by its zeroth component we would obtain: 

A + t dh(El 

AH = *fc • (2.36) 

Of course, this is not simply the complex conjugate of f)2.34p . as in the BPS case W > 0, whilst 
for the non-BPS case W < (recall h = y/fw\). 



2.2 Solutions to the Attractor Equation 

It is relatively straightforward to obtain solutions to the above attractor equation when we 
restrict ourselves to a system with only DO and D4 branes (go and p a charges). This attractor 
was found in [5] by minimizing Vbh - in fact, the solution here is not restricted to STU, i.e. 
there may be any number of the moduli. We consider here, though, the complementary solution 
with DQ and D2 branes (p° and q a charges). The superpotential is then given by: 

W(z\ z 2 , z 3 ) = q a z a + p°z 1 z 2 z 3 . (2.37) 

It is also straightforward to write down the symplectic invariant I\. 

h = V±4p°gig 2 g 3 • (2.38) 

A positive sign under the square root corresponds to the non-BPS attractor and a negative one 
to the BPS one, since, as discussed above, W(r) is generally positive for charges corresponding 
to BPS solutions and negative for those corresponding to non-BPS ones (it is evident from the 
form of the solutions below that p ^^^ < for BPS attractors p gi<?2<?3 > otherwise). 

^From (12.381) it is trivial to use (I2.36P to obtain the BPS attractive moduli values as: 



1 ■ / 92^3 2 ■ / 91 93 3 • / 9l?2 (n qm 

u p°qi V V <72 V V Q3 



3 We are grateful to E. Gimon who suggested to look for a more general form of eq. (2.34) which would 
include (D0-D6) system. 
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The non-BPS (from ( 12.341) ) results have a similar form: 



1 • Q2qs 2 • QiQa 3 ■ QiQ2 (n An , 

V p u qi V V Q2 V P 93 

These results have been confirmed by minimizing Vbj? and by using the form of the attractor 
equation given in [TB] - these details of those calculations can be found in appendices A.l and 
A. 2. It's worth noting at this juncture that care must be taken with these solutions to ensure 
that their signs are such that e K is positive. This point is discussed in somewhat more detail in 
the appendix. We also observe that (I2.34p and (12.361) can be used to find the attractor values 
of the moduli for the D0-D4 system addressed in [5] and [IE]. As expected the BPS answers 
are: 



i • IQoP 1 2 • <loP 2 3 ■ doP 3 fn A ^ 

p z p 6 y p^p 1 y p L p z 



And the non-BPS: 



¥s • - 2 = -h - q 4l ■ = "M ■ (2-42) 

p l p 6 y p 6 p L y p L p z 



As mentioned above, for a non-BPS attractor we are required to confirm that the extremal 
point of the potential does correspond to a minima. It turns out that this is a somewhat subtle 
and involved calculation. The details can be found in appendix A. 3, but the short-answer to 
the question of stability is (as in the D0-D4 case discussed in [5]): "It's stable." 

Before moving on to a discussion of the complete solutions for the moduli fields we summarize 
our knowledge of STU black hole attractors: 

• All extremal STU black holes can exhibit attractor behavior, with the values of the moduli 
at the horizon found by extremizing the effective potential Vbh- 

• Only those systems with extremum of the potential a minimum will form stable attrac- 
tors. These systems include those with Y = (q ,0,0,p a ) and those with T = (0, g a ,p°, 0). 
Unstable "attractors" include those with all four types of charge, T = (qo,q a ,P°,P a )] this 
result was established in [3] (see appendix A. 3 for details). 

We shall return to the issue of stability when we consider multi-centered non-BPS black holes 
in section 4. Now, however, we move on to the general solution for the moduli. 



3 The Exact Non-BPS Solution 



We now demonstrate explicitly the construction of exact solutions for the fields in the non-BPS 
attractor. We begin with a general discussion, and then write down and prove the form of the 
solutions. 
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3.1 Harmonic Functions and U-Duality 



So far, we have calculated the value of moduli at the horizon of the supersymmetric and non- 
super symmetric STU black holes via the attractor mechanism. It turns out, however, that we 
can do better than that - we can find solutions to the equations of motion for the moduli, 
allowing us to obtain their values everywhere. This has already been done for the BPS case 
(details can be found in [17], [18], [19] and [27]). and we will proceed in an analogous fashion. 

We work with the single-center spherically-symmetric, static metric ansatz: 

ds 2 = -e 2U dt 2 + e~ 2U dx l dx l . (3.1) 

The basic idea is that one takes the horizon (r = r^) values of the moduli, z a (q A ,p A ) and 
replaces the charges with harmonic functions: 

H A A dh{H) 

- A (^W) = ^—mSk ■ ^ 



H°-i 



dH 

The metric is found through: 

e- 2U = h(H) (3.3) 

The harmonic functions H are: 

H(r) = (H A , H A ) = (r A , r A ) r + (h A , h A ) , r = . (3.4) 

Ii is the symplectic invariant defined above and there is a constraint on h: (h, T) = 0. H arises 
from the solution to the vector equations of motion where the vector fields can be represented 
by the doublet (F A ,G A ). These fields are not independent: 

G A = Re Af AE (z, z)F A - ImAT AS (z, z) * F A (3.5) 

The equations of motion then require a vector field potential given by (^4 A , B A ) where F A = dA A 
and G A = dB A . The vector fields are related to the harmonic functions by: 

^mn = ^mnpdpH , G Amn = —€ mn pdpH A . (3-6) 

The timelike components of F and G are fixed in terms of the spatial ones through ( 13. 5j) . 

That this result holds in the BPS case has been established explicitly; there is, however, no 
such proof for non-BPS black holes. Instead, we shall apply the above algorithm to generate 
an ansatz which can be checked explicitly using the STU equations of motion. Before we do 
so, though, we offer an argument as to why one might expect the above prescription to work. 

There exists a manifest [SL(2, Z)] 3 symmetry of the equations of motion which mixes the 
Maxwell equations and Bianchi identities: 

A A (x)\ _ ( A B\ ( A A (x)\ 

B A (x) ) \C D { B A (x) ) [6J) 
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Here ( ^ ) is a global Sp(8, Z) matrix such that A T C - C T A = 5 T L> - £> T £ = and 



n C D 

A T D — C T B = 1. Only a subgroup of this symmetry, corresponding to [SL(2, Z)} 3 duality 
symmetry, is unbroken at the level of solutions to the equations of motion. The details of 
the embedding of [SX(2,Z)] 3 into Sp(8,Z) are given in [21] and [26]. Simultaneously with the 
duality transformations of the vector potential doublet the 3 moduli also undergo a standard 
fractional transformation of the type: 

W*>Y = aJ ^^ (3.8) 
cz\x) + a 

This transformation of the moduli follows from the consistency of equation ( 13.5f) and leaves the 
metric invariant 

(9iw(x))' = 9iw(x) , (3.9) 

allowing us to generate from one solution all the others related through duality. Clearly, such a 
symmetry places strong constraints on the form of the moduli fields. Further, our explorations 
of functional forms have been unable to find any solutions with the appropriate symmetry other 
than those given by (13.21) . when we simultaneously impose the obvious criteria that the moduli 
give the correct horizon values, remain bounded everywhere and have derivatives that do the 
same while also vanishing at the horizon (see [16] for the origins of these restrictions). 



3.2 The Proof 

We shall begin by establishing that the appropriate ansatz solves the D2-DQ brane system and 
then find the general (almost - see below for details) solution through a symmetry argument. 



3.2.1 The D2-DQ System 

The four dimensional stationary, spherically-symmetric effective Lagrangian of the Maxwell- 
Einstein action (derived from (12.21) ) is: 

C(U(r),z a (r),z a (r)) = (u 2 + G al z a i h + e 2U V BH ) . (3.10) 

r is the inverse of the radial coordinate (such that the horizon is at r = — oo and the derivatives 
denote differentiation with respect to r. The gravitational and scalar equations of motion 
derived from the above Lagrangian (13.101) are 

U = e 2U V BH , (3.11) 
z a + T a bc z b z c = e 2U G al d- b V BH . (3.12) 

In addition, there is also a constraint on the system: 

U 2 + G al z a i l - e 2U V BH = c 2 , (3.13) 
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where c 2 = 2ST = for extremal black holes. 



As described above, to generate solutions to these equations we take the horizon values of 
the moduli, (12.401) . and replace the charges with appropriate harmonic functions: 



-2U 



2 V / # #1#2#3 



H 2 H 3 
H°Hi 
H\H 3 

H\H 2 



—i 



—i 



(3.14) 



Substitution followed by explicit calculation of the terms in the gravitational equation of motion 
gives: 



Ue- 2U = V BH 



%-y/H°H 2 Hz 



4 



(#i) 3 / 2 



4 



(H 3 )^ 



(# 2 ) 3 / 2 

(p°) 2 

(#0)3/2 1 



H°H!H 3 



H\H 2 H 3 



(3.15) 



In a somewhat more involved calculation we can also verify that the ansatz solves the scalar 
equations of motion. Picking the 1 direction, we have: 



^ + ri 1 (i 1 ) 2 



e w G ll diV BH 

e 2U+K G n (g 22 {D- x D 2 W)D 2 W + G 3 ~ 3 (DiD 3 W)D 3 W + 2{DiW)W )(3.16) 



Once again substitution and tedium give the result we desire: 
z' + T^z 1 ) 2 = e w G ll d l V BH 



i fH 2 H 3 ^- 3/2 



2 



qi(H»H 2 H 3 ) 2 - {jP)\HiH 2 H.. 



Finally, we check the constraint: 



U 2 + G al z a f = e 2U V BH 



"I" T-I-O T-I-O I 



HI HI H 2 (IP) 



(3.17) 



(3.18) 



3.2.2 The Generic System: (D0,D2,D4,D6) 



To generalize the previous result, we need to establish that the prescription of the previous 
section, i.e. the substitution of harmonic functions for the electric and magnetic charges in the 
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values of moduli at the horizon, solves the equations of motion when the moduli and metric 
are given by the expressions above: 



z A (H(x)) 



H 



A _ - dh{H) 
1 dH A 



-2U 



H° 



; dh(H) 
' dH 



(3.19) 



Given the hugely increased complexity of I± when all charges are present (not to mention Vbh) 
an explicit calculation would painful. Accordingly, we're going to do something else. 

The symplectic invariance of special geometry ensures that the Lagrangian of our theory 
has an Sp(8, Z) symmetry, which reduces to [SL(2, Z)] 3 at the level of the equations of motion. 
This symmetry group can be used to generate solutions for generic charges by rotating the 
expressions for the moduli in the D2-D6 system (eq. (13.141) ) obtained above. To do this we 
take the following element of [SL(2, Z)] 3 : 



a b 
c d 



1 
1 



1 
1 



(3.20) 



in which ad — be = 1. Under this transformation, the charges and harmonic functions of the 
D2-D6 system transform in the following way 



f H°\ 

H 1 
H 2 
H 3 
Ho 
Hi 
H 2 

Further, we know that under transformation (13.201) . the moduli coordinates transform as 

- 1 ' b , 



p 1 
p 2 
p 3 



/ dp \ 

bp 

cq 3 
c<? 2 
-bqi 
dq x 
aq 2 
\ aqs ) 



( dH° \ 
bH° 
cH 3 
cH 2 
-bHt 
dHi 
aH 2 
\ aH 3 ) 



(3.21) 



z 1 



az 



z 2 



z 3 



, , - - ~' 5 - (3-22) 

cz 1 + a 

The solution of equations of motion for the case of generic charged can be swiftly obtained 
from fl3T22|) . 



To check whether (I3.19P (with a set of generic charges) coincides with (13.221) it is sufficient 
to establish that (13.221) will solve the equations of motion. To do this we first we observe, as 
proved in [17], that the hyper-determinant W is invariant under (I3.20p . In our particular case 
this is simple to see from: 



W(H) i-> W{H) = -4(ad - bcfH^HxH^ = W(H) 



(3.23) 



4 We observe that the charges (and relevant harmonic function) are not completely independent and we 
have these two relations among them: p 3 ^ = p 2 q2 and p^qi — —p°qo. If we also make nontrivial SL(2,Z) 



transformations on the second and third moduli this dependence disappears. On the other hand, the above 
analysis shows that this procedure can be applied for successive transformations. Therefore, we do not lose the 
generality of the argument. 
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Next, using (I3.19p . we obtain the modulus along the first direction of the moduli space for the 
generic case as 



^ _ H l v^W| + i{H\H 2 H 2 + H 3 H 3 - H°H - H 1 ^) - 2H°H 2 H 3 ) 



(3.24) 



H°^/\W\ + i{2H l H 2 H 3 - H°(H A H A )) 



This can be reexpressed in terms of the harmonic functions of the D2-DQ system: 



by/H°H 1 H 2 H 3 - iaH 2 H 3 



a[—i 




az 1 + b 



(3.25) 



d^/WHJhH 3 - icH 2 H 3 



ci—i 




cz l + d 



where we have used the fact ad— be = 1 repeatedly. This is in complete agreement with f)3.22p . A 
similar check confirms that z 2 = z 2 and z 3 = z 3 , as required. Hence, we can conclude that (13.191) 
gives the complete solutions to the equations of motion for stationary, spherically-symmetric 
non-BPS attractors. There is one small caveat, what we have actually done is produce a general 
solution from a restricted subset of D2-D6 attractors with ho = h 1 = h 2 = h 3 = 0; while we 
believe that this simplification will not effect out result (it amounts to adding a constant 
electric/magnetic potential, something which clearly does not affect the vector equations of 
motion), we have been unable to check the equations of motion for the unconstrained D2-D6 
case. 



So far we've managed to establish that many of the properties of STU non-BPS attractors are 
analogous to those of their BPS counterparts. In particular, along with the basic attractor 
behavior established in [I], we have also demonstrated explicitly (at least for the STU model) 
that complete expressions for the moduli can be formed from the attractor values, by replacing 
charges with harmonic functions. Given this success, it seems opportune to push a little further. 

It has been shown in [ITj-fTH] and [2Z|-[22! th & t; static, spherically-symmetric black holes are 
not the only BPS attractors one can find. Rather the attractor mechanism can be expanded to 
stationary spacetimes (dropping the requirement of spherical symmetry) with the metric: 



The most general situation described by this geometry is one consisting of multi-centered, 
charged black holes. Each black hole sits at a point x a and carries a set of charges T a . The 
vector fields are then fully described by a set of harmonic functions (a fact that doesn't rely on 
the BPS properties of the black holes, or the number of centers): 



4 Split Non-BPS Attractors 



ds 2 = - e 2l/(x) (dt + uM) 2 + e- 2U{x) dx i dx i 



(4.1) 




(4.2) 
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Bates and Denef have demonstrated (in [27]) that for the BPS case these harmonic functions 
can be used to construct exact solutions in an analogous fashion to the single center black holes. 
Ui is defined through the harmonic functions: 

Ui = (H, dH) . (4.3) 

As before there is an integrability condition: 



1=1 



The above fixes the distance \xj — Xj| unless the charges are mutually local, with (1^,1^) = 
- this is the case when Ui = 0. We also have a similar constraint on h as for the single center 
attractor: (h, T) = 0, where T = J2i IV m this setup the BPS solutions are given by: 



Z (#(*)) , . dh[H ) 



i- 



dH 

-2U 



e~ w = h(H). (4.5) 
Where we define the symplectic invariant I\ as before: 

h{H) = \Z(H)\ 2 + \DZ(H)\ 2 = V\W(H)\ > , (4.6) 



For the non-BPS STU model we thus conjecture that, again analogous to the single center 
black hole, solutions are given by: 

,r : (//(.r)) - 



H°-i 



dh(H) 



dH 

e~ 2U = h(H) (4.7) 

Of course, here W(H) < 0. Our principle motivation here, once again, is the underlying 
duality symmetry and the constraining restrictions it must place on the form of any solution. 
In theory, we can check this solution in an entirely similar fashion to the spherically-symmetric 
cases discussed above. However, the process is somewhat more involved, and as yet it has not 
proved possible to carry out the required calculations. 

There are, though, good reasons to believe that multi-centered, non-BPS, extremal black 
hole attractors exist; and that the behavior of the moduli fields is described by (14.71) . As 
discussed in [4 J and [TH] the attractiveness of black holes is a result of their extremality. Specif- 
ically, extremal black holes have a near horizon geometry described by a Bertotti-Robinson 
(AdS2 x S 2 ) product space with an infinite throat. It is this infinite throat that leads to attrac- 
tor behavior and it should thus be clear that if extremal multi-centered stationary solutions 
can be found, then each black hole will have the same near-horizon behavior as it would do 
in isolation. Accordingly, we should expect each horizon to be an attractor. Furthermore, it 
is apparent that solutions for z a of the form (14.71) will reduce to the single-center ones in the 
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near-horizon limit - suggesting that they may well be the appropriate global expressions for the 
moduli. The missing piece in this argument, and the reason one would want to check equations 
of motion explicitly, is that we cannot be certain that stationary, multi-centered solutions exist 
for extremal, non-BPS black holes. 

Before moving on to a discussion of our results, there is one last observation that we should 
make. Denef has observed in [28] that for a multi-center (charges T a ) BPS black hole the 
behavior far from all centers should be that of a spherically symmetric solution with charge 
T = ^ a r a , leading to a split attractor flow. The same should hold true for non-BPS case, 
though here we could have an interesting additional effect. Whilst BPS extremal black holes 
are always stable, non-BPS ones are not - so although systems with Ti = (q Q , 0, 0,p a ) and those 
with T 2 = (0, q a ,P°, 0) give stable attractors, systems with T = Ti + T 2 = (go, <la,P ,P a ) do not. 
Thus we might expect that apparently unstable "attractors" correspond to stable multi-centered 
configurations when examined at appropriately short distances from the horizorjE 



5 Discussion 



In the preceding sections we have discussed some of issues arising in non-BPS STU attractors. 
In particular we have found that the system consisting of p° and q a charges can give rise to a 
stable non-BPS attractor - a result that was obtained using both the attractor equation (in 
two different forms) and through the minimization of an effective black hole potential. Further 
we have demonstrated that one can construct exact expressions for the moduli, the metric and 
the vector fields using harmonic functions through a prescription analogous to that used for 
BPS black holes. 

Additionally, we have argued for the existence of stationary non-BPS attractors, with an 
associated split attractor flow. This allows us to make the following conjecture regarding 
unstable STU non-BPS "attractors": 

Conjecture. The apparent unstable nature of non-BPS STU "attractors" can be explained 
through the existence of a stable multi-center solution that is only resolved at sufficiently short 
distances. It is only far away from both centers that the gradient flow of the system appears 
repulsive. 

We hope to present a proof of this conjecture in future work. 

Finally, we note that while everything we have discussed only applies to STU black holes, 
we see no reason in principle why the behavior we have found and the arguments we have made 
for both exact solutions and multi-center attractors will not apply to more general extremal 
non-BPS black holes. 

5 Thcrc is some additional subtlety here - since the gradient flow will have to switch direction when the flow 
splits we should expect the line of marginal stability ([29], [2H]), where the split takes place, to have somewhat 
different properties to that in the BPS case. 
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A Further Details of the D2 — D6 system 

In this appendix we discuss various details of the D2-DQ system - developing the solution 
established above by confirming it through some alternate methods and analyzing its stability. 



A.l Minimizing the Effective Potential 

Although we calculated the values of moduli at the black hole horizon by using the attractor 
equation we can, instead, directly extremize the effective potential; i.e. solve d a VBH = 0: 

d a V BH = e K (G b£ (D a D b W)D 5 W + 2{D a W)w) . (A.l) 

Since the superpotential only includes odd powers of the coordinates for the system of D2-D6 
brane, we assume that the ansatz of the minimization of the black hole potential equation is 
purely imaginary z a = —z a (see [1]). Now, we want to compute the ingredients of (lA.ip . The 
first covariant derivative of the superpotential is: 



W(z\ 



z 



3> 



D a W = (d a + d a K)W = Z J_ fa — • ( A - 2 ) 



Note that in the above expression, z a is substituted in the argument of the superpotential for 
z a . With the ansatz z a = —z a , (1A.2I) then gives: 

DiW = wiW) = w-2 (q ^ + q ^) (A 3) 

DlW = wklllfl = YLl^k+^fl , (A.4) 

z 2 — z 2 2z 2 

W(z\z 2 ^) _ W - 2(q 1 z 1 + q 2 z 2 ) 
DsW ~ z*-z* 2? ' (A ' 5) 

When we drop the arguments of the superpotential they should be understood to be z 1 , z 2 and 
z 3 . Before we consider the non-SUSY solutions, we find the supersymmetric ones, D a W = 0: 



' p u qi V P 92 V P 1 3 



and also 



z =-iJ — 5— ,z =-iJ — 5— ,z = -U — 5—- A.7 
V P <?i V P I2 V P <b 



As briefly discussed above, we should examine these solutions to see whether they respect the 
positivity of the Kahler potential - a necessary condition to ensure positive kinetic terms in the 
action. A simple check shows that the first SUSY solution flA.6j) violates the positivity of the 
Kahler potential and therefore, it is ruled out. Thus only the latter (IA.7I) is acceptable. We 
should also notice that these solutions will only exist when p°qiq 2 q 3 < 0. 

The next object to calculate is the second covariant derivative of W: 

D a D b W = (d a + d a K)D b W - T c ab D c W . (A.8) 

For our particular STU model: 

D a D b W = (1 - 8*) W ^'''jfl*'" ) . (A.9) 

(z a — z a )(z b — Z b J 

Notice that z a and z b substitute z a and z b in the argument of W . For the z a = —z a ansatz of 
our D2-DQ brane system, the above reduces to: 



W — 2(n z a 4- rihZ b ) 

D a D b W = (1 - 5 ab )- 4 ° a 6 ( no sum on a and b ) ■ ( A - 10 ) 



Now, we are ready to form the equation ofan explicit expression for the minimization of the 
potential d a VBH = 0. Using (12. 8p . flA.3j) . and flA.101) . we find the following set of equations: 

(W - 2{q l z 1 + q 2 z 2 )^j (w - 2(q lZ l + g 3 ^ 3 )) +w(w~ 2(q 2 z 2 + g 3 ^ 3 )) = , (A.ll) 
(W - 2(q 1 z 1 + q 2 z 2 )^ (w - 2{q 2 z 2 + q 3 z 3 )^ + w(w - 2{q l z l + q 3 z 3 )^ = , (A.12) 
(W - 2(q 1 z 1 + q 3 z 3 )^ (w - 2{q 2 z 2 + q 3 z 3 )^ + w(w - 2(q 1 z 1 + q 2 z 2 )^ = . (A.13) 

It is clear that fl A. 7[) satisfy these equations. Now, we want to find the non-SUSY solutions 
(D a W 7^ 0). If we divide the above expressions by one another, we obtain: 



(W - 2(q 1 z 1 + q 2 z 2 )) 2 = (w - 2(q 2 z 2 + q 3 z 3 )) 2 = [w - 2(q 1 z 1 + q^ 3 )^ . (A.14) 
The complete solutions to these relations are: 



V p"qi V p ?2 V p <?3 

If we consider all these possibilities for (IA.14j) . then we find eight sets, which can be categorized 
into two groups: 

A : {(+, +, +), (+, -, -), (-, +, -), (-, -,+)}, (A.16) 
B : {(-, -, -), (-, +, +), (+, -, +), (+, +, -)} . (A.17) 

In above, each parenthesis shows the sign of (z 1 , z 2 , z 3 ) respectively. It turns out that the 
elements of group A violate the positivity of e K and therefore this group is unacceptable. 
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However, no such problem exists for group B, and therefore its elements provide acceptable 
non-SUSY solutions. Furthermore, the last three elements of group B are physically equivalent 
solutions in which the sign of two fields are the same, but opposite to the third. It's worth 
noting that that the attractor equation gave us only the (— , — , — ) solution; however, this point 
has little effect on the rest of out argument. 

Finally, we note that for all non-SUSY solutions, we must have p ^^^ > 0. This in turn 
means that we cannot simultaneously have both SUSY and non-SUSY solutions. 



A. 2 An Alternate Attractor Equation 



In this section, we find the values of the moduli at the horizon of the black hole by solving the 
attractor equation in the form given in [16]: 



T = 21m 



2Z 



(A.18) 



As before T is the set of magnetic and electric charges T = (p A , g^) and II is the covariantly 
holomorphic period vector: 



IT = e K l 2 



( 1 \ 

\F a J 



(A.19) 



We have already calculated most of these ingredients, with our only remaining task to calculate 
the covariant derivative of II: 



/ 



DAI = e K ' 2 



d d K 



\ 



5 a d + (d d K)z a 
d d F + (d d K)F 
V 9 d F a + (d d K)F a ) 



(A.20) 



It is straightforward to calculate the two terms of the attractor equation for our ansatz z a 
These are: 



2Im(Zn) 



W 



4z 1 z 2 z 3 



( 1 \ 






z 2 z 3 
z l z 3 
\ z l z 2 ) 



21m 



{D- a D- b Z)G ac G M D c Z 
2Z 



D d Ji 



Az 1 z 2 z 3 W 



( Y \ 





Y lZ 2 z 3 
Y 2 z 1 z 3 
\ Y^z 2 ) 



(A.21) 
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We have defined Y , Yi, Y 2 , and Y 3 to be: 

Y = W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) + W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) + W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) 
Y l = W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) + W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) - W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) 
Y 2 = W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) + W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) - W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) 
Y 3 = W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) + W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) - W(z\ z 2 , z 3 )W(z\ z 2 , z 3 ) . 

(A.22) 

Forming the attractor equation flA. 18[) . we find the following: 

W 2 - Y = Ap°Wz 1 z 2 z 3 (A.23) 
W 2 + Y a = 4q a Wz a (no sum on a) . (A.24) 

As is clear, we have four equations, but only three of them are independent. This can be easily 
seen if we add the four to one another and use F x + Y 2 + Y 3 — Y Q . To solve these equations we 
first note that Yi — Aq x Wz l = Y 2 — Aq 2 Wz 2 = Y 3 — 4q 3 Wz 3 and observe this leads to: 

qtz 1 = ±q 2 z 2 = ±q 3 z 3 . (A.25) 

If we take qiz 1 = q 2 z 2 = q 3 z 3 we find: 

1 P°3l/ 1N 3 A ( J. , Al/ 1N3 



z y _ r*L( z y ^ + ^(z 1 ) 3 = . (A.26) 
Q2q3 J V J 

We immediately realize that the first parenthesis produces the supersymmetric solution: 



z 1 = ±iJ-^ , (A.27) 
V P Qi 

The second parenthesis of (IA.26I) then gives us the non-supersymmetric attractor: 



zi = ±ij^. (A.28) 



As we see, different choices of ()A.25j) produce all the different SUSY and non-SUSY solutions 
that we found in the previous section. Recalling that the positivity of e K requires we choose a 
specific set of signs for the above solutions leads to (1A.7I) for the supersymmetric solution and 
elements of group B (see IA. 16|) for the non-supersymmetric ones. 



A.3 Stability 

In the previous section, we found both SUSY and non-SUSY attractors. The former are always 
stable, since the second derivative of the black hole potential is proportional to the Kahler 
metric. However, for non-SUSY solutions, we need to check stability explicitly. In order to do 
this, we need to find the full mass matrix. 
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First, note that at the extremum of the potential we have d a d b V B H = D a D b V BH ; therefore, 
we can use the covariant rather than ordinary derivative. Then the holomorphic-holomorphic 
and holomorphic-antiholomorphic parts of the mass matrix are given by 



D a D b V BH 
D- a D b V BH 



e K ( G cd (D a D b D c W)D 3 W + 3(D a D b W)W ) 



icd 



e K [ - n a b - ac G cE (D d W)D s W + G cd (D b D c W) (D- a D d W) 



all 



+G- ab G cd (D c W)D d W + 3(D- a W)D b W + 2G- ab WW ) , 



(A.29) 



(A.30) 



where TZ bac is the Riemann curvature tensor of the moduli space. We begin with the holomorphic- 
holomorphic piece. This requires us to compute the third covariant derivative of the superpo- 
tential: 



D a D b D c W = -\e abc \ 



W(z\z 2 ,z 3 ) 



(z a - z a )(z b - z b )(z c - z 
Using this result and defining B ab = D a D b V BH \ horizon we have: 



, (no sum on a, b, and c) . (A. 31) 



B = lp°y/iP 



5i525a 




(A.32) 



To find the antiholomorphic-holomorphic piece of the mass matrix, we need the Riemann cur- 
vature of the moduli space: 



TZ a 



no sum on a 



(A.33) 



(z a - z a ) 2 ' 

All other components of the Riemann curvature tensor vanish. Taking this result we can crawl 
through a long computation to obtain A ab = D a D b V B H as: 



horizon 



/ 2qi_ J_ J_ \ 



5i5253 



<?2 93 93 
1 2q 2 



92 

1 



V 



93 9193 91 
1 1 2g 3 



/ 



92 91 9192 

Hence, the full mass matrix of the black hole potential is given by: 



(A.34) 



M 



A B 
B A 



/ g9i J_ 

292 



P°5i5253 



9293 

1 



93 

1 



9193 
1 



92 


J_ 

93 

V - - 

\ 92 91 



93 





92 

J_ 

91 

2g 3 



9192 
1 



92 
J_ 

91 






j_ 

93 
1 



2gi 



9293 
1 



93 

1 



92 



93 



1 



91 
1 



292 

9193 
J_ 

91 



i ^ 

91 





92 
J_ 

9192 



/ 



(A.35) 



Note that this will give rise to a term in the effective potential which takes the following form: 

z ] Mz where, z = (z 1 , z 2 , z 3 , z 1 , z 2 , z 3 ) . (A.36) 
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To explore the stability of the non-SUSY solutions we need to diagonalize the mass matrix. 
This is straightforward and the eigenvalues are given by: 

nn /~6 P° qi /~6 /~~o p ° q3 /~o f p ° qi i p ° q2 i P° q A\( a ? 7 \ 

0, 0, VP u ?i?2«3 , VP 9i?293 , V P Q1Q2Q3 , VP 9i9293 1 1 MA.37) 

9293 9i93 9i92 v 9293 9i93 9i92^ J 

There are two flat directions, with the rest of the eigenvalues positive. For the former we find 
that the associated eigenvectors are: 

-^,0,1,-^,0,1] and (--,1, 0,-— ,1,0] . (A.38) 
9i 9i / V 9i 9i / 

Both of these (as is clear from (IA.36I) ) are in a real direction - moving along potential in the 
direction of the eigenvectors sends z a — > z a + Sx a , where 5x a is real. 

In order to explore the stability of our extrema with respect to the flat directions, we need to 
examine higher order derivatives. Recall that we initially proceeded with covariant derivatives 
instead of ordinary derivatives (since DDVbh = cPVbh at extremal points of Vbh)- But this 
is no longer true at higher orders, so one needs to revert to the flat derivative. Unfortunately, 
this becomes rather unpleasant, as we forced to break the covariant form, hence we choose to 
proceed in a calculationally more tractable fashion: We expand the black hole potential around 
the non-supersymmetric solutions and then concentrate on the behavior of cubic and quartic 
terms of the expansion. If the coefficients of the cubic terms are non-vanishing, then we can 
conclude that the non-supersymmetric solution is unstable. However, if the coefficients vanish, 
then we need to consider the quartic terms. If these terms are positive for flat directions, 
independent of the values of the parameter of the expansion, then the non-supersymmetric 
solution is indeed stable. 

First, using (12. 8p and (IA.30 . we notice that we can rewrite the black hole potential in the 
following way 

Vbh = e K (\W(z 1 , z 2 , z 3 )\ 2 + \W(z\ z\ z 3 )\ 2 + \W(z\ z 2 , z 3 )\ 2 + \W(z\ z 2 , z 3 )\ 2 ^j . (A.39) 



Now, we want to expand the above potential around the non-supersymmetric solution (IA.15I) 

as 



1 ' ,12 ■ /9l93 .2 3 • /9l92 |3 f K AC\\ 

V V 9i V P 92 V P 93 

where e a is a real parameter of the expansion. Why only real? Well we are, of course, only 
concerned with the quartic expansion in the flat directions, and we have already established that, 
for our first non-supersymmetric solution, these are real. Accordingly Vbh has the following 
expansion 

V BH (z a ) = V BH {z%) + 0((e a ) 2 ) + 0((e a ) 3 ) + 0{{t a f) + • • • , (A.41) 

where Zq is the value of moduli for non-supersymmetric solution. Note that the flat directions 
in (1A.38P correspond to the plane q^e 1 + q 2 e 2 + q 3 e 3 = in notation of (1A.40D . 
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The information in 0((e a ) 2 ) is encoded in the mass matrix which has already been calculated 
and we need to know about higher order terms along flat directions. Explicit computation shows 
that the cubic terms vanis 10 C((e a ) 3 ) = 0. Thus we must examine the quartic terms and we 
find: 

0((e a Y) = 4e^ 2 oyf^(e 1 6 2 ) 2 + ^(6 2 e 3 ) 2 + ^( e 1 e 3 ) 2 

^ 93 Ql <?2 

+ gi (6 1 ) 2 e 2 e 3 + g2 e 1 (6 2 )V + g3 6 1 6 2 (e 3 ) 2 ) 
= 2e K ^^^[a 2 1 + a 2 2 + al + (a 1 + a 2 + a 3 ) 2 } . (A.42) 

Here a\ = qi^qzC a>2 = qi^q^ and 03 = gie 1 ^ 2 '■ As we see, no matter how the parameters 
of expansion change, 0((e a ) 4 ) is manifestly positive. 

This implies that the non-SUSY solution (1A.15|) is indeed stable. Finally, we mention that it 
can be (and has been) checked that the other non-SUSY solution (— , +, +) also gives a positive 
mass matrix. 
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